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ABSTRACT 


The GOES satellites launched in the 1980s are carrying an instrument called the VISSR At- 
mospheric Sounder (VAS), which is designed to provide temperature and moisture profile sound- 
ing capability foi mesoscale weather systems. As a prelaunch study of this capability, VAS radi- 
ance fields are simulated for pre-thunderstorm environments in Oklahoma to demonstrate three 
points: (1) significant moisture g.'fl jients can be seen directly in images of the VAS channels, 

(2) temperature and moisture profiles can be retrieved from VAS radiances with sufficient accu- 
racy to be useful for mesoscale analysis of a severe storm environment, and (3) the quality of 
VAS mesoscale soundings improves with conditioning by local weather statistics. 

Even thougli the simulated VAS soundings have the usual limitations in absolute accuracy 
and vertical resolution (especially in the lower tropospheric moisture retrievals), it is still possible 
to discern potentially unstable mesoscale structures from the synthetic VAS radiances for a partly 
cloudy pre-thunderstorm environment. The RMS tropospheric profile errors are ±1° C and ±25 
percent in temperature and mixing ratio, respectively. The results represent the optimum retriev- 
ability of mesoscale information from VAS radiances without the use of ancillary data. The sim- 
ulations suggest tliat VAS data will yield the best soundings when a human being classifies the 
scene, picks relatively clear areas for retrieval, and applies a “local” statistical data-base to resolve 
the ambiguities of satellite observations in favor of the most probable atmospheric structure. 
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VISSR ATMOSPHERIC SOUNDER (VAS) SIMULATION EXPERIMENT 
FOR A SEVERE STORM ENVIRONMENT 

1. INTRODUCTION 

Present atmospheric soundings from operational satellites are designed for the collection of 
temperature and humidity data for global analysis and for input to numerical models that simulate 
the large-scale general circulation. The soundings are performed from sun-synchronous polar or- 
biting satellites, which provide complete coverage of the globe with a spatial scale of several hun- 
dred kilometers twice per day per satellite. The entire system, which includes the sounders on 
the satellites, the algorithm?: for reducing the data and retrieving the desired atmospheric variables, 
and the general circulation models, is aimed at improving the intermediate-to long-range forecasts 
of the large-scale circulation patterns. 

The need for understanding short-lived weather phenomena, such as severe local storms, and 
tlie growing requirements for more accurate local weather forecasts that require atmospheric data 
with a finer spatial scale and more frequent observations point toward a different satellite system: 
one centered around soundings from a geosynchronous platform. Severe local storm research and 
prediction require frequent observations at subsynoptic to mesoscale resolution in order to cap- 
ture the relevant changes in the atmospheric mass and wind fields. The incorporation of ancillary 
data (e.g,, conventional surface and radiosonde reports) from the meteorologically active areas is 
expected to become an important aid to the satellite sounding process. Likewise, an understand- 
ing of the physics of the phenomenon under observation should be helpful in modeling the radi- 
ation transfer for the retrieval algorithms. 

The VISSR Atmospheric Sounder (VAS), which is now operating as an integral part of the 
Geostationary Operational Environmental Satellite (GOES) system, is designed to meet the require- 
ments of observing subsynoptic to mesoscale temperature and moisture distributions. Ultimately, 
the goal is to observe, understand and predict the development of severe storms. The geosynchro- 
nous VAS will watch tlie development of mesoscale weather systems, typically surveying the 
United States with 12 spectral bands once per hour at 15 km (nadir view) resolution. Satellite 
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soundings derived from such multispectral scenes vy'ill have 30 to 90 km effective resolution, since 
several pixels must be used to average radiometric noise and to deal with broken cloud cover in 
the net sounding field of view (SFOV), 

At this time, the polar orbiting satellite sounders are operationally useful over the data-sparse 
oceans where both the underlying uniformity of the ocean surface and the passive microwave data 
arc important data processing aids for synoptic scale retrievals. These soundings show some sys- 
tematic error patterns, such as larger errors near the surface and tropopause, biased performance 
in partly cloudy areas, and low sensitivity to horizontal gradients (Phillips, et, al,, 1979). 

Operational satellite soundings over land have similar errors: the subsynoptic gradients are 
underestimated, and there are significant errors where the underlying topography is unusual or 
persistently cloudy (Schlatter, 1981). Operational satellite soundings are usually ignored in favor 
of the conventional radiosonde network, despite the inadequate space-time resolution of the con- 
ventional synoptic data for mesoscale developments. Case studies over mesoscale networks with 
satellite radiances from the Vertical Temperature Profile Radiometer (VTPR) have been encourag- 
ing. VTPR retrievals at 70 km resolution over the Atmospheik Variability Experiment (AVE) 
produce temperature fields and precipitable water estimates comparable to the AVE data (Hillger 
and Vender Haar, 1 977). VTPR radiance correlations over the National Severe Storms Laboratory 
(NSSL) network show spatial structure comparable to the meteorological correlations in the NSSL 
analyses, although the retrieved moisture information is not well determined (Hillger and Vonder 
Haar, 1979). 

Polar orbiting sounder studies are hampered by the lack of unambiguous, independent, colo- 
cated “ground truth” data, taken with space-time resolution as good as the satellite observations. 
This VAS simulation study avoids these problems by declaring that selected mesoscale profile anal- 
yses are “ground truth,” and then retrieving profiles from corresponding synthetic radiance fields 
so that unambiguous estimates of retrieval accuracy can be made, The quantitative results provide 
a benchmark for mesoscale retrieval accuracy available from the VAS sounder for optimal sound- 
ing conditions without help from ancillary data or from physically modeled retrieval algorithms. 
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Tliis paper (lescril)cs a simulation experiment designed to estimate the accuracy of VAS 
soundings within a severe weatlicr environment, given cloud*free observing conditions. The pre- 
launcli experiment is specifically designed to determine: 

1 . Wiiether moisture gradients can be seen directly with the VAS images, 

2. Whether temperature and moisture profiles and gradients can be retrieved from VAS radi- 
ances with sufficient accuracy to be useful for mesoscale analysis of a severe storm en- 
vironment, 

3. Whether the quality of VAS soundings improve v/ith conditioning by local weather sta- 
tistics. 

A brief summary of the VAS instrument and the VAS channels is provided in section 2. The 
experiment is defined in section 3. Simulated imagery results are discussed in section 4 with a 
statistical analysis of the soundings presented in section 5, The application of VAS soundings in 
the analysis of precursor conditions for severe convective storms Is presented in section d, and a 
discussion of the results follows in section 1 . 

2. DESCRIPTION OF VAS 

During the 1970s, the Visible Infrared Spin Scan Radiometer (VISSR) was carried on the 
SMS/GOES series of geosynchronous satellites. VISSR infrared images in the 1 1 micron window 
are formed by stepping the scan line each spin by its own width (7 km resolution at nadir). Dur- 
ing the 1980s, the GOES satellites will carry an upgraded instrument called the VISSR Atmospher- 
ic Sounder (VAS). The VAS instrument lias programmable features which allow a ground-based 
operator to interact with the satellite by using the most recent observations to command when, 
where and how VAS will acquire data. The user will be limited only by the available time, filter/ 
detector combinations and radiometric constraints (see Table 1 at the end of this technical mem- 
orandum). The geosynchronous location can be exploited to observe the development of meso- 
scale weather systems. For example, VAS can take sounding quality (low noise) radiances at 15 
km (nadir view) resolution along a full earth east-west swath with a north-to-south extent large 


enougli to cover Oklahoma every 20 minutes or with a north-to*south extent large enougli to 
cover tlie conterminous United States every hour. 

2.1 The VAS fnstnmmit 

VAS lias 12 tlicrmal infrared channels between 4 and 15 microns. The channels are chosen 
to distinguish the effects of tropospheric temperature, moisture and cloud cover upon the upwell- 
ing radiances. All 12 channels can be operated at 15 km resolution (instantaneous geometrical 
field of view at nadir), and channels 4, 5, 7, 8, 9 and 10 can be operated at 7 km resolution. 

(The other channels are too noise- or diffraction-limited to be operated at higlier resolution.) The 
VISSR racUometer has been upgraded to yield higli precision values (10 bits) from six detectors, 
carefully calibrated with outer space (3® K) and an internal hot target (320® K). Radiometric 
noise is reduced by “dwell averaging” several spins on the same line and channel; the number of 
spins allocated to the channels is termed the “spin budget,” Noise can also be reduced after data 
collection by judicious averaging of the clear adjacent fields of view, 

Table 1 summarizes key design features of tlie 12 VAS channels: central wavelength, wave- 
number and bandpass; purpose for sounding; main absorbing gasses and other important influences; 
single-sample noise; radiometric noise requirements for sounding; and a nominal spin budget for 
dwell averaging each channel. Lower spin budgets can be used to trade radiometric accuracy for 
speed and latitude coverage. In that case, pixel averaging of the 15 km (nadir view) instantaneous 
geometrical fields of view (IGFOVs) into, say, a 60 km sounding field of view (SFOV) is possible 
to regain radiometric accuracy at the expense of detailed resolution. Nominal calibration speci- 
fications are ±1.5® K absolute and ±0,5® K relative brightness temperature AT* in all channels 
for a typical earth scene. For aomparison, the previous VISSR instrument had calibration errors 
of rouglily ±3® K. 

2.2 VAS Channel Sensitivities 

VAS transmission functions, r (n, p), are computed for each channel (designated with wave- 
number v) at 40 pressure levels p, from ,2 to 1000 mb. The transmissions are done for COj, Nj, 


NjO, NO, CII 4 , Oj, plus HjO, using. Hnc-by-line absorpucm coefficients (McClatchcy* et. a!., 

1973) and important molecular continuum absorption cross sections (Burch, et, a/., 1970*1; Rob* 
erts, ct. (tin 1976). Corresponding weightlnp functions, d r/dln P, are shown for the U,S, Stan- 
dard Atmosphere In Figure 1, which Is presented at the end of this technical memorandum, for 
the temperature sounding channels and in Figure 2 for the moisture and window channels, Tlic 
short wave channels are indicated with dashed lines. Tlie redundancy between the long and short 
wave channels will be used to detect cloud contamination by exploiting the differential sensitivity 
of the Planck function to cl.'vud-top temperature. The sharp peaks of tlm moisture channels in 
Figure 2 arc related to tlie exponential average distribution of water vapor. Consequently, the 
altitude of tlicir peak response is not very predictable, since the radiance accumulates from only 
the top few millimeters of tlie total water burden. 

Figures 1 and 2 show tlnvr llie lower troposphere (500 to 1000 mb) is not well rcsolvctl by 
any one of the VAS channels. The radiation from the temperature and moisture structures just 
above the boundary layer (700 to 900 mb) is mixed with both surface and mid-tropospheric radi- 
ances in all the key VAS channels (5 to 8 and 12). The “man in the loop” concept will be an 
important factor to successful sounding with VAS (Smitli, at. al, 1978) where radiances will be 
confused by low clouds, topograpliy and irregular surface emission. Also, tlie inclusion of surface 
temperature and moisture information into the retrieval sclieme can be of great value in properly 
determining tlie atmospheric boundary conditions (Hayden, 1980). This prelaunch sounding ex- 
periment will test whetlier the lower tropospheric VAS cliamiels can independently detect potenti- 
ally unstable temperature and moisture structures using simulated clear column radiances over a 
featureless level surface which is thermally coupled to tiie bottom of the troposphere. 

Table 2 summarizes some calculated responses of tlie VAS cliannels to the clear atmosphere: 
the pressure level of the peak of each standard weigliting function; the pressure range over whicli 
the middle 80 percent of the radiance accumulates; tlie net sensitivity of the brightness tempera- 
ture to a r C increase in the entire air temperature profile or in the dewpoint profile or in the 
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surface temperature alone; mean and standard deviations of VAS simulated brightness tempera- 
tures for a set of 32 radfiCf^ondes distributed over Nortli America; and “signarVnoise estimates 
made by comparing these variances to the single-sample noise listed in Table 1. Roughly speak- 
ing, the ‘‘signai’" (variance) in the tropospheric channels is 5 percent of the total brightness, and 
the noise Is 0.5 percent giving a “slgnar/noise of about 10 to 1. 

Tlie net sensitivity estimates in Table 2 suggest simple channel combinations for special en- 
vironments. One promising combination is VAS channels 7 and 8, which form a “split window” 
for water vapor dctermiriation In the lower troposphere, since the two channels have different 
sensitivity to moisture variations at nearly the same wavelength (Prabhakara, at, at., 1974). For 
instance, a SFOV with 5” C higlier dewpoint than tlie U.S. Standard will appear 1.8“ K colder at 
12 microns due to the increased absorption by water vapor, but only 1.0“ K colder at 11 microns. 
Consequently, this “split window” should be able to directly display horizontal moisture grad- 
ients, such as pre-thundsrstorm “dryline*’ environments with warm diy air next to cool moist air. 

3. experiment design 

The VAS simulation experiment is designed to determine the potential accuracy of VAS 
temperature and moisture retrievals within a realistic severe weather environment and the sensitiv- 
ity of these retrievals to localized statistical conditioning. 

Since VAS will be directed to take radiance observations over a developing mesoscale event 
(e.g., thunderstorm or huiTicane), the VAS sounding algorithms should be prepared to take advan- 
tage of the information implicit in the classification of the scene by the VAS user. This study 
demonstrates the retrieval advantage enjoyed by a user who has prepared for pre-thunderstorm 
soundings by the use of objective statistics gathered from previous mesoscale radiosonde analyses. 
While it is difficult to quantify the accuracy of a human classification, one can measure the post- 
classification retrieval advantage by comparing soundings made with appropriate “local” training 
to soundings made with nnselective “global” training. 



3.1 NSSL Database 

Temperature mid moisture profiles for eight severe local storm events were selected from tiie 
radiosonde network operated by the National Severe Storm Laboratory (NSSL) in 1976 (see Ta- 
bles 3 and 4) to be utilized for VAS radiance simulations. Tiie storm environments range from 
mild to chaotic gradients, most having “dryllne” characteristics with weak temperature gradients 
but strong moisture gradients (Mostek, et. at., 1980). Ten different “frames” (times) were anal- 
yzed from the eight events* so that each frame could be used to simulate corresponding VAS radi- 
ance observations. Tlie events were generally analyzed early in the storm development, wliich 
are relatively cloud-free, representing times when the VAS should be able to operate effectively 
as a sounder. VAS’ nominal resolution (15 km at nadir view) will be roughly 20 km at the lati- 
tude of Oklahoma, 

Tiie NSSL radiosondes (launched every 90 minutes from sites rouglily 50 km apart) were in- 
terpolated onto a 21 x i 1 grid with approximately 20 km resolution using a Barnes objective 
analysis scheme (Barnee, er, al„ 1971). Thegriddcd profiles serve as a “ground truth” for the 
sounding simulation experiment. Figure 3 is an SMS 1 1 micron image of the Texas-Oklahoma 
area at the time of the 2030 GMT 29 May 1976 frame listed in Table 4. The analysis grid is 
overlaid on Figure 3, showing a 200 x 400 km area encompassing the NSSL network. In order 
to simulate the coarser SFOV coverage (50 km) expected from VAS under operational conditions, 
a 7 X 5 subset of the 21 x 11 grid is used for clear air radiance simulations of all 10 NSSL frames, 
This 7x5 subgrid (50 km resolution) is the set of SFOVs actually used in the statistical and grad- 
ient studies in the following sections, while the original 21x11 grid (20 km resolution) is only 
used for image simulation. The SMS image in Figure 3 shows us that we can expect large enough 
cloud-free areas on either side of a developing instability for VAS to carry out soundings. 

3.2 Statistical Design 

Tlie NSSL profiles and the corresponding simulated radiance fields are separated into two 
sets, dependent and independent, for establishing and testing the correlations between atmospheric 
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structure and VAS radfances. 1'lic first 6 of the 10 frames are used as the dependent training set, 
determining the clear alr-to-radianee correlations needed for constructing the first guess and regres- 
sion matrices used In tliis study. The remaining 4 of the 10 NSSL frames (drawn from just 2 
NSSL events) are used as the independent test set for calculating the retrievability of atmospheric 
information from clear SFOV radiances, Since each frame contains 7x5 SfiOVs, the NSSl.6 
training set has 210 dependent points, and the NSSL4 test set has 140 independent points, for 
retrieving temperature atui moisture at 1 5 tropospheric levels, Both subsets have similar mean 
and covariance statistics. 


The quality of profile retrievals is mcasiued by the RMS residuals between tlie simulated 
VAS soundings and tlie “ground mitli” profiles drawn from the NSSL analyses. Tiie quality of 
gradient retrievals is judged from horizontal and vertical contour plots of analyzed soundings com- 
i)arcd to the “ground trutli” NSSL analyses. Tlie usefulness of simulated VAS soundings for meso- 
scale analysis is demonstrated by the retrieval of a convectively unstable air mass. 

The impact of progressively more detailed statistical conditioning is studied by comparing 
RMS profile residuals from retrievals conditioned with two other dependent sets; a “global” set 
supplied by NOAA from 32 widely distributed radiosondes and a wet/dry subset of the 6 “local” 
dependent NSSL frames. The NOAA32 arc drawn iTom North American radiosonde/rocketsonde 
sites with a bias towards unusual profiles (Smith and Woolf, 1976). The “global” set is used lierc 
to cstablisli statistics tliat are so broadly selected In time and space as to be unclassified. The 
NSSL wet/dry sets were selected by sub-classifying tlie individual NSSL6 profiles as “wet” or 
“dry,” clepcnding on whether the mixing ratio at 850 mb was more or less tliang gm kg"’, respec- 
tively. 


4. VAS RADIANCE SIMULATION 

In order to simulate a VAS frame over Okiahoiiui, radiances are computed from transmission 
functions for the temperature and moisture profile at each NSSL SFOV, 
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4J Radiation Transfer 

The VAS radiance simulations use line-by4ine transmission functions, r p), computed for 
each temperature, T(p), and mixing ratio, Q(p), profile, Radiances R(p) in each channel, at nom- 
inal wavenumber v, are computed by radiation transfer of the Planck function, B(j^, T), up from 
a perfectly emitting surface; 

T (p) = r {v, p; T(p), Q(p); lines, continuum), (1) 

R = B (T(s))t(s) + /° B (T(p))rir(p), 

s 

The U.S. Standard Atmosphere is used to fill in the stratoplieric profile above the last re- 
ported radiosonde level, Topographic irregularities are accounted for by setting the 1 000 mb and 
surface temperatures equal to the first reported radiosonde level, and surfaces are 100 percent e- 
missive, giving isothermal bottomside structures with no radiance contrast. This assumption re- 
sults in radiances tliat arc correct for the actual topography, but it also results in a strong correla- 
tion between the lower troposphere and the surface, so that tlie lower tropospheric retrievals rely 
very lieavily upon the window radiances, 

Cloudy radiances are modeled with a simple cloud fraction n(x, y) at pressure level p ' with 
cloud-top temperature T(c) equal to the air temperature T(p') for the SFOV: 

R = I ' - n] [B (T (s)) T (s) + / B (T (p)) dr (p) ] (2) 

s 

+ iiB (T (c)) r (c) B (T (p) ) dr (p). 

c 

4.2 Simuiated VAS Images 

Figure 3 shows the actual SMS 1 1 micron picture of the Texas-Oklahoma area taken at 2030 
GMT 29 May 1 976, In order to make realistic cloud height and coverage estimates, the cloud 
brightness temperatures T* (VISSR) are matched to a temperature profile estimate T^(p) of the 
scene, with iteration to produce a field of variable cloud fraction n(x, y) with the same cloud top 
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pressure p', Tliese cloud estimates are combined with the transmission runctions for the analyzed 
SFOVs in the 21x11 grid to produce cloudy radiance fields, which are converted into VAS grey- 
level images with 20 km resolution. Some of these images are shown in Figures 4A, B and C 
where significant moisture gradients arc visible through the openings in the clouds. 

Figure 4 shows three simulated images of partly cloudy channels: the “VISSR window” 
(channel 8 at 11.2 microns), the “wet window” (channel 7 at 12.7 microns), and the difference 
in their brightness temperatures (T*(Ch8) - T’*'(Ch7)), The cloud tops appear as dark, cold 
SFOVs in the radiance images. In the cloudless areas, the brighter areas are warmer and dryer 
and the darker areas are cooler and wetter due to greater water vapor absorption. The appearance 
of this relative moisture gradient is enhanced by taking the difference between these “split win- 
dow” clianncls. The differential image in Figure 4C has the brightest areas in the western por- 
tion of the NSSL frame where the two-channel differences arc tlic greatest due to the absorption 
by the excess jnoisturc tlicrc. The maximum contrast shown in Figure 4 is 5.5“ K, which should 
dominate the expected single-sample noise (NEAT = ±0.3“ K in Table 1). Figures 1 1 and 12 
show the actual cross sections of the temperature and moisture fields retrieved from this frame. 
These are discussed at length in section 6. 


4.3 Correlatioti Between Moisture and VAS Radiances 

The moisture gradient visible in the simulated VAS images suggests that we should look at 
the correlation between the VAS “split window” radiances and the moisture content of the lower 
troposphere. A clear correlation would be a basis for subdividing the local NSSL profiles into 
wet/dry subclasses for use in finely tuned dryline retrievals. 

Figure 5 plots the brightness temperature differences (“ K) between the “split window” 
channels against the mixing ratio (gm cm“^) at 850 mb for the six dependent NSSL cases. There 
is the expected trend in the “split window” differences, with wetter SFOVs being darker due to 
the moisture blocking upwctling radiance more effectively at 12 microns. However, the scatter 
from fratnc to frame blurs tlic simple difference to moisture correlation enougli to foil its use ns 
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a direct predictor of absolute moisture content. A plot of total precipitable water (gm cm“^) 
against the briglitness difference (® K) in the “split window” produces roughly the same scatter 
aa seen in Figure 5. The scatter in these simulated radiances is due entirely to natural variations 
in the temperature and moisture fields at other levels. Radiometric noise would contribute only 
another ±0.3® K to the brightness temperature differences. Correlation plots like Figure 5 under- 
state VAS* net sensitivity to moisture, since each channel is influenced simultaneously by several 
atmospheric parameters. Retrieval algorithms combining all 1 2 VAS channels should be used to 
separate variables and to improve vertical resolution. 

5, RESIDUAL ERRORS IN SIMULATED VAS SOUNDINGS 

Tile conventional figure of merit (FOM) for profile retrievals is the RMS residual profile be- 
tween the sounding and “ground truth.” The “ground truth” in this study is simply the tempera- 
ture and moisture profiles from which VAS radiances are simulated. 

5.1 Conditioned Regression Retrievals 

A minimum-?".formation, linear regression algorithm is used to retrieve temperature and mix- 
ing ratio profiles from the simulated VAS radiances. This technique (Fleming and Smith, 1971) 
is chosen for the prelaunch demonstration of the VAS capabilities because it has enough flexibil- 
ity to adjust both to VAS’ programmable noise levels and to the statistical conditioning available 
in historical mesoscale radiosonde collections. Such statistical conditioning should improve re- 
trievals by resolving ambiguities of vertical structure in favor of the most probable structure ex- 
pected from objective experience (Chesters, 1980). 

Regression retrievals assume that the radiance response 6R (v) in a channel at wavenumber 
n, is a linear combination of atmospheric variations, 6X(p), where X(p) is temperature, moisture, 
or other state parameter in the atmosphere as a function of pressure, p. Consequently, an esti- 
mate of the atmospheric variation, 5X'(p), can be made from a corresponding linear combination 
of variations in the observed radiances: 


II 



6R 00 « A (V, P) • 5X (p). 
5X' (p) « C (p, »0 • 5R (p). 


( 3 ) 


Tlie symbolizes matrix iivaUlplication with terms being summed over matching indepemlent 
variable values. With many observations of X and R, a retrieval matrix, C(p, p). can be estimated 
froni the radiance to atmosphere cross correlations in a set of observations: 

C (p. P) « < 5X (p) 5R (P') > • I< 5R (p') 5R (p) > . (4) 


The symbolizes a statistical average over a sample learning set of simultaneous, colocated 

observations, and the “f,,.)’*® ” symbolizes matrix inversion, The basis of the <6X5R> correlation 
is the presence of T and Q within r and B in equation I. Likewise, the <5R6R> channel-to- 
channel cross correlation is based both upon tlic underlying correlations within the state variables 
X ami upon the overlap among the channels dr/dlnP (see Figures I and 2). 

Direct application of the simple regression algorithm in equation 4 to a small set of available 
satellite and radiosonde ilata usually leads to retrieval matrices with large cocflicients and oppos- 
ing signs. This over-interprets the small lluctuations in real data. The mininium-intbrmation 
scheme conditions retrievals against this effect by incori)orating an estimate of the radiometer 
noise in each channel <e (p)*>. 

Mathematically, the atmospheric variations and the corresponding radiance changes for a de- 
pendent set of events arc linearized about the mean, and a single-sample noise, e(p), is included: 


X = <X> + 5X. 

R = <R> + SR + e, 

5X' = C • [5R -f e] , 
d<[6X' - aX]2>/dC = 0, 

C = <5X6R> • K6R5R> + <e2>]“' 

Equation 5 uses least squares techniques to solve for a matrix C (p, v), whicli julnimizes tlie un- 
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certulnty in the retrievals 6X' (p) averaged over a training set of dependent events. The radio* 
meter noise can be expressed directly in radiance units (erg cm“- sec*"* wavenumber"* steradian**, 
abbreviated “erg/etc,*’), or it can be computed from the single-channel variance <5R (p)^> divided 
by a corrcspojiding “signar’/uoise ratio 7 (p): 

<e2> = <6R2>/72 (6) 

Note that the “signal" is the variance in the radiance, not the total brightness of a channel. The 
retrievals in this study are all done with a single “conditioning factor," 7 = 100 / 1 , which is opti- 
mistic for most of the VAS channels (sec Table 2), 


5.2 "Locor vs. 'Vlobar Retrieval Quality 

Simulated radiances for tlie dear SFOVs are used to determine the accuracy of temperature 
and moisture retrievals from soundings roughly 50 km apart, A “local" retrieval matrix is calcu- 
lated using the six NSSL sets of 7 x 5 SFOVs by correlating tlie radiances with tropospheric tem- 
peratures and mixing ratios. For comparison, anotiier retrieval matrix is calculated from the 
NOAA32 set of widely distributed radiosondes and their simulated radiances. Retrievals are done 
on the four independent NSSL test frames of 7 x 5 SFOVs with both the locally and widely con- 
ditioned matrices. 

The RMS residuals in Figures 6 and 7 arc the information not retrieved from these 140 clear, 
independent SFOVs because of vertical smoothing and noise conditioning. The NSSL retrievals 
arc Indeed much worse when they are made with matrices trained by statistics from a set of 
widely distributed radiosondes, The residual errors in the "globally" conditioned NOAA32 retri- 
evals of the independent NSSL SFOVs are about ±4° C in temperature and ±25 percent in mois- 
ture, which are quite inaccurate for mesoscale weather analysis. The residual uncertainties in the 
“locally" conditioned NSSL 6 retrievals are about ± 1 “ C for the tropospheric temperaUire profile 
and ±25 percent for the mixing ratio profile, which are well within the VAS sounding require- 
ments, although the moisture uncertainty is still too large to be useful for absolute mesoscale 
analysis. However, Figure 13, discussed in section 6 , shows that local conditioning retrieves more 
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accurate relative moisture gradients, which arc not discernible in these composite RMS profile 
statistics. 

The moisture residuals are poorly determined in the lower troposphere, since VAS has 
only one moisture sensitive channel there, as indicated by the weighting functions in Figure 2, 
Neither the 100 mb nor the 1000 mb retrievals are meaningful, since the original NSSL radio- 
sondes were rarely reported above 200 mb or below 920 mb. Profile uncertainties arc dominated 
by mean residuals (not shown) over several unresolved vertical layers, as is normal with passive 
infrared sounding instruments. The temperature residuals are normally larger at the tropopause 
because there is little radiance iVom these cold, unresolved layers. Both the moisture and tem- 
perature residuals are normally greater near the surface because of the naturally greater atmospher- 
ic variability there. 

Under real operating conditions, the uncertainty in the lower tropospheric profile will be 
even greater because of the unresolved clouds. The window channels will detect uneven emission 
from a surface that is less ccjrclatcd with the “surface” air temperature than is assumed in this 
simulation. Operational temperature errors are expected to be about t2‘’ C; the temperature re- 
siduals in this study are only ±1“ C because the radiances are simulated for optimal conditions, 

5.3 Fifie Tuning the "LocaT' Retrieval Statistics 

The better quality of the “locally” conditioned soundings encouraged us to try even more 
detailed conditioning, “Dryline” conditioned retrievals are attempted by computing additional 
regression matrices from “wet” and “dry” subsets of the six dependent NSSL frames with SFOVs 
classified by whether the low level moisture, Q(850 mb), is more or less than 8 gm kg"‘ . Retrie- 
vals are done upon selected SFOVs from the four independent NSSL test frames. The relatively 
bright (warm and dry) SFOVs in the 1 1 micron window are retrieved with the “dry” matrix, and 
the relatively dark (cool and moist) SFOVs are retrieved with the “wet” matrix. The RMS pro- 
file residuals from these wet/dry conditioned retrievals are plotted in Figures 8 and 9. For com- 
parison, the original NSSL6 “locally” conditioned residuals are also plotted in these figures. The 
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wet/clry training scarcely improves the residuals, indicating that we have reached the point of di- 
minishing returns by sucli detailed local conditioning. 

5.4 Retrieval of Vertically Integrated Parameters 

Since ’passive satellite radiances in each channel are accumulated over several hundred milli- 
bars of the troposphere (see Figures 1 and 2 and the third column of Table 2), one can expect 
satellite retrievals to represent vertically integrated parameters more accurately tiian single-level 
parameters. We examine the residuals in two integrated parameters derived from retrieved pro- 
files for tile four independent NSSL frames: 

1. Precipitable water, measured in gm cm“^', calculated by integrating the mixing ratio re- 
trievals from 100 to 920 mb. 

2, Geopotential thickness, measured in meters, calculated by integrating the temperature 
retrievals from 500 to 920 mb. 

Table 5 lists both the mean and the standard deviation in the “ground truth” set and the 
residual errors for regression retrievals done with the “global” NOAA32 and with the “local” 
NSSL6 training sets. Note that the RMS error is the combination of the mean and standard de- 
viation in the residuals and that a successful retrieval must have RMS errors less than the variance 
(standard deviation) in the “ground truth” set. Table 5 shows that: 

1. The “ground truth” precipitable water in the four NSSL test frames is 2.49 ±0.39 gm 
cm"^. The “global” NOAA32 retrievals did not improve the information about the pre- 
cipitable water, since the RMS error (±0.39 gm kg”* ) is as large as the original variance. 
The “local” NSSL6 retrievals improved the uncertainty in the total precipitable water 
from ±0.39 to ±0,19 gm cm"^, or only an 8 percent error in the total. The mean errors 
in the precipitable water retrievals are reasonably small for both retrieval matrices, which 
indicates tliat the “split window” (VAS channels 7 and 8) is working as planned. 

2. The “ground truth” geopotential thickness in the four NSSL test frames is 4998 ±28.5 m. 
The NOAA32 “global” retrievals actually degraded the information about the thickness 
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field, with an RMS residual of ±31,2 m, due mostly to a mean underestimate of -30.5 m. 
The errors in the “globally’’ retrieved thickness field reflect the corresponding underesti- 
mates in the “globally” conditioned temperature profile retrievals discussed in the prev- 
ious section. The NSSL6 “locai” retrievals greatly improved the information about the 
thickness fields, leaving an RMS residual of only ±3,3 m. 

In summary, the residual errors for the four independent NSSL test fields indicate that total pre- 
cipitable water should be retrieved with, less than ±25 percent error (even with poor statistical 
conditioning) and that gcopotential thickness fields should be retrieved with “locally” conditioned 
methods in <».der to avoid serious mean errors. Analyzed fields of retrieved precipitablc water 
and geopotential thickness will be examined in the following section. 

6. MESOSCALE ANALYSIS OF VAS SIMULATED SOUNDINGS 

An important objective of tills VAS simulation experiment is to determine if V.4.S tempera- 
ture and moisture retrievals can adequately delineate a convectively unstable air mass prior to 
the outbreak of severe convective storms. In order to demonstrate this, “ground truth” fields 
are compared to retrieved vertical cross sections, horizontal cross sections, and contour maps of 
vertically integrated parameters. 

The 29 May 1976 case is selected for this part of the study, since this severe weather cut- 
break is an excellent example of a “dryline” propagating into Oklahoma. The lack of cloud cover 
over the NSSL network at the chosen time provided an opportunity for the realistic use of clear 
column radiances without cloud corrections, simulating optimal sounding conditions. The Na- 
tional Weather Service (NWS) surface analysis at 0000 GMT May 1 976, shown in Figure 1 0, de- 
picts a moderately strong low pressure system centered over the Texas panhandle. The most no- 
table feature near the NSSL network is the dewpoint gradient due to the drier air moving east- 
ward toward southwest Oklahoma. The SMS infrared image for 2030 GMT 29 May 1976 in Fig- 
ure 3 shows a relatively clear re^on within the special network ahead of the developing convec- 
tive stonn cells in the Texas panhandle. 
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6.1 Analysis of Horizontal Cross Sections 

Tlic ’’ground truth” 850 mb temperature, mixing ratio, and equivalent potential temperature 
(fljj) fields at 850 mb and 2030 GMT 29 May 1976 are shown in Figures 1 lA, D and G, This 
case lacks a significant temperature gradient across the network with only a 3° C range within 
the domain. On the other hand, there is a moisture maximum in western Oklahoma within which 
explosive convective development occurred by 0000 GMT 30 May 1976. The 0^ analysis reflects 
this moisture maximum and the resulting gradients across the central portion of the network. 

Figure 1 1 also shows the 850 mb temperature, mixing ratio and 0^ analyses derived from the 
VAS simulated retrievals using the 7x5 SFOVs at 50 km resolution. Contour plots are shown 
for the two sets of retrievals, conditioned by the NOAA32 dataset and the NSSL6 dependent 
dataset. Comparison among the analyzed fields in Figure 1 1 reveals a general similarity in the 
main features with significantly different details; 

1. The retrieved temperatures are within ±2° C of the ’’ground truth,” However, the 
NOAA32 retrievals have a spurious temperature minimum in the center of the network, 
reversing the actual small east-west gradient in the western portion of the domain. The 
NSSL6 retrievals more accurately depict the structure of the warmer air, although the 
magnitude of the gradients is still weaker than the “ground truth,” 

2. The retrieved 850 mb mixing ratios are significantly underestimated by both sets of 
soundings shown in Figure 11. However, the NSSL-conditioned retrievals yielded mixing 
ratios 2 gm kg“Miigher than the NOAA32 retrievals in the region of maximum moisture 
and better accounted for the gradients of the moisture field, with the gradients in the 
central portion of the network nearly three times larger. There is a corresponding im- 
pact upon the gradients in the analyses. 

In summary, Figure 1 1 shows that the “locally” conditioned retrievals do a better job of retriev- 
ing the horizontal structure in both temperature and moisture. The location and relative structure 
but not absolute magnitude, of the low level 0^, gradient are well retrieved. 
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6,2 Analysis of Vertical Cross Sections 

Figure 12 shows corresponding vertical cross sections of temperature, mixing ratio and 0^ 
constructed from west to east througli the middle of the network, intersecting the region of max- 
imum moisture and convective instability (as measured by dO^/dz). The cross sections of the 
“ground truth” data, in Figures 12A, D and G, again reveal the small temperature gradients and 
the maximum in moisture and 0^ in the western portion of the networks Of particular interest 
is the location of the strong vertical gradient of the mixing ratio just above the 850 mb level 
and the corresponding vertical gradient in 0^. The large negative value of (\0Jc\z in the western 
portion of the domain illustrates that this region is indeed primed for the severe convective storms 
that developed there within the next 3 hours. Examination of the vertical cross sections of the 
NOAA32 and NSSL6 retrieval analyses reveals these interesting differences: 

1. Both retrievals underestimate the moisture content and the vertical gradients. The 
NOAA32 cross section in Figure I2E actually has a larger moisture value in tlie lower 
troposphere, but shows a more gradual vertical gradient than the NSSL6 cross section in 
Figure 1 2F. However, the NSSL6 moisture retrievals better depict the magnitude and 
location of tlie vertical moisture gradient. Tins behavior can be related to tlie makeup 
of the moisture profiles in the NOAA32 and NSSL6 training datasets. The lower tropo- 
spiieric moisture information embodied in tlie NOAA32 is determined by tropical pro- 
files, which have a large concentration of moisture in tlie boundary layer and a gradual 
decrease with height. The moisture information in the NSSL dataset is basically drier, 
but typical profiles have a sharp gradient near the 850 mb level. Consequently, these 
pre-thunderstorm structures are only approximately reproduced by the corresponding 
retrieval matrices, because the broad weighting functions shown in Figures 1 and 2 are 
unable to resolve the vertical structure of the moisture in the lower troposphere. 

2. Both retrieval sets in Figure 12 underestimate 0^ in the lower troposphere due to the 
deficiencies in the moisture retrievals. Yet, both retrievals do in fact capture a region of 
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inaxlnuim 0^ and dd^/dz In the region of niaxlinum convective instability, Tiic NSSL6 
retrievals locate the maximum dd^/dz gradient more accurately at the critical 850 mb 
level tlian do the NOAA32 retrievals, Tims, the "locally” retrieved 0^ field is more rep- 
resentative of tlic actual convective instability sliown in Figure I2G. 

3. The lifted index (LI), which is a measure of the atmosplieric instability, is computed for 
the most unstable region of the NSSL network, ns determined by the 0^ cross sections in 
Figure 12G, The LI Is -10 for the "ground trutii” data, indicating a very unstable en- 
vironment conducive to the development of severe convective storms. For botli tlie 
"global” NOAA32 retrievals and the "local” NSSL6 retrievals, the LI is -4. Tims, the 
simulated VAS retrievals capture a potential instability, although its magnitude is reduced. 
In summary, Figure 12 shows tiiat the "local” conditioning does a better job of estimating the 
vertical moisture structure. Tlic potentially unstable air mass is well located in the retrievals, but 
its magnitude is underestimated due to the lack of vertical resolution in the VAS ciiannels and to 
the lack of ancillary surface data in this simulation. 


6.3 Analysis of Integrated Parameters 

The RMS statistics for prccipitable water and geopotential thickness retrievals over all four 
independent NSSL test sets indicate that these integrated parameters should be a better measure 
of satellite sounding quality tlian single-level cross sections. To demonstrate this point for an 
analyzed field, contour maps of total water and thickness arc produced from the retrievals for 
the 2030 GMT 29 May 1976 frame. Figures 13 A, B and C sliow prccipitable water (100 to 920 
mb), and Figures 13D, E and F show gcopotential tliickness (500 to 920 mb). Comparison of 
the "ground truth” fields to the "global” NOAA32 retrievals and the “local” NSSL6 retrievals 
confirms the trends seen in t!ie 850 mb level cross sections in Figure 11: 

1. The “ground truth” precipitable water field in Figure 13A has a strong gradient from 
east to west (ranging from 1.47 to 3.25 gm cm“^ with a mean and standard deviation 
of 2.27 ±0.46 gm cin"^). Tlie precipitable water is somewhat underestimated with both 
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retrievals, althougli the residual error is small (averaging +0.09 ±0.37 gm cnr^ for tlie 
“global” field, and only -0.03 ±0.12 gm cm“^ for the “local” field). The location and 
gradients in the total water arc roughly captured in both retrievals, but with noticeably 
better detail in the “local” retrievals. Apparently, the “split window” channels can de- 
liver enough unambiguous information about the total water content to avoid serious 
retrieval errors even with poor statistical conditioning. 

2. The “ground truth” gcopotcntial thickness field in Figure 13D has only weak gradients 
from east to the southwest corner (ranging from 5018 to 5048 m with a mean and stan- 
dard deviation of 5029 ±5.2 m). The small variance reflects the weak underlying temper- 
ature gradients. The “globally” conditioned thickness retrieval field in Figure 13E is a 
very poor representation of the “ground trutli’h the weak eastern minimum is exagger- 
ated, the spurious warm retrievals in the center produce a spurious maximum, the gradi- 
ent in the soutliwest corner is missed, and all thickness values are much too low (aver- 
aging -33.5 ±5.9 m). Apparently, the radiance signature in the moisture channels is 
being incorrectly interpreted as tropospheric temperature variations by the “global” 
regression statistics, since Figure 13E looks like a combination of the moisture and tem- 
perature features seen in both Figures 13A and D. The “localiy” conditioned thickness 
field in Figure 13F is in much better agreement with the “ground truth” field, with most 
of the error in the form of a small overestimate for the entire “locally” retrieved field 
(averaging only +4.0 ±2.8 m), 

In summary. Figure 13 shows that analyzed fields of parameters integrated over the simulated 
VAS retrievals are estimated well enough to distinguish outstanding mesoscale features, especially 
when the retrievals are “locally” conditioned. Residuals are dominated by mean errors in the 
absolute strength of the field. Therefore, VAS can be expected to be a very good instrument 
when used for capturing relative horizontal gradients within a mesoscale field. 
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7. SUMMARY AND DISCUSSION 


The VISSR Atmospheric Sounder (VAS) is designed to operate as an integral part of the 
GOGS system, providing an interactive sounding capability from geosynclironous orbit for the 
first time, VAS lias 12 calibrated tliermal infrared channels in the COj and ll^O molecular bands 
between 4 and 15 microns. Tiie channels arc cliosen to distinguish the effects of tropospheric 
temperature, moisture and cloud cover upon the upwelling radiances. The geosynchronous nature 
of VAS is a key clement for tlie observation of the tlcvclopment of viiesoscalc weather systems, 
typically surveying the United States in all 12 bands once per hour with 15 km (nadir view) res- 
olution. Satellite soundings derived from such multispeetral scenes will have 30 to <50 km effec- 
tive resolution, since several pixels will have to be averaged to reduce ridiometric noise and to 
deal with broken cloud cover in the net sounding field of view. 

Tills report describes a prelauncii sounding simulation experiment to demonstrate: 

1. The visibility of lower tropospheric moisture gradients in the VAS radiances. 

2. The detectability of a potential convective instability in a pre-thunderstorm environment, 
using optimal soundings from the VAS channels, 

3. The sensitivity of residual errors in VAS regression retrievals to “local" statistical condi- 
tioning for the kind of mesoscale phenomenon under investigation. 

In order to estimate the accura of tiie three-dimensional VAS soundings within severe local 
storm systems, 12-cliannel radiance fields arc simulated for thunderstorm environments observed 
in Oklahoma. Temperature and moisture profiles for 10 time frames arc drawn from the radio- 
sonde network operated by the NSSU in 1 976. The profiles are gridded using a Barnes analysis, 
covering Oklahoma with a grid at 20 km horizontal resolution and the troposphere at 1 km verti- 
cal resolution. Radiances are simulated for SFOVs in each frame, using line-by-line transmission 
functions for the temperature and moisture profiles. 

Simulated VAS images of partly cloudy pre-storm conditions are shown for one of the NSSL 
frames wliere a corresponding SMS 1 1 micron picture is available. Radiances are computed at 
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full VAS resolution (20 km projected view per SFOV). Grey level images in the “split window” 
cliannels at 11 and 12 microns show that realistic cloud cover does not seriously mask the under- 
lying precursor storm conditions, Moisture gradients are visible in tlic clear areas witli contrast 
greater than tlie expected VAS noise levels. 

Temperature and moisture profile retrievals arc done at the expected SFOV resolution by 
using a coarser (50 km) SFOV subset drawn from the finer (20 km) grid, These simulated radi- 
ances arc separated into two sets; a “local” dependent training set used to determine clear air- 
to-radiance correlations and an independent test set used to determine the retrievability of atmos- 
pheric information from clear SFOVs, The impact of detailed statistical conditioning upon a 
regression retrieval scheme is demonstrated by comparing these “locally” conditioned retrievals 
to soundings done with two other dependent training sets: a “global” set derived from 32 North 
American radiosondes and a wet/dry subset of the “local” dependent NSSL frames. For the lat- 
ter comparison, the fields are partitioned into a “wet set” and a “dry set,” using a threshold mix- 
ing ratio of 8 gm kg“* at 850 mb. 

Results from the simulated retrievals include: 

1 . Temperature profile residual errors in the troposphere are ± 2 ° C across a 3° C gradient, 
and corresponding mixing ratio residuals are ±2 gm kg’*' across an 8 gm kg“' gradient. 

2. Horizontal and vertical temperature and moisture structures arc retrieved in the lower 
troposphere, and potentially unstable air is recovered from simulated VAS radiances hav- 
ing the resolution and noise limits expected from the instrument. The absolute magni- 
tude of the retrieved potential instability is significantly underestimated by the remote 
sounding process, which reflects the lack of vertical resolution in the 12 VAS channels. 
Relative temperature and moisture patterns are retrieved, although the magnitudes of the 
gradients are likewise underestimated. 

3. Retrievals done witli regression matrices conditioned with “local” weather statistics have 
smaller residual errors than those conditioned with “global” data. These results are due 
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to the better resolution of ambiguities in vertical structure in favor of the most probable 
structure derived from a “local'’ statistical database. A furtiier subdivision of the depen- 
dent NSSL6 dataset into “wet” and “dry” subsets did not make any statistically signifi- 
cant improvement in the “local” retrievals, 

4. VAS should be capable of retrieving vertically integrated meteorological parameters. The 
main moisture and airmass features are captured in the form of retrieved fields of prccipi- 
table water and geopotential thickness. These simulations had only ±0.2 gm cm“^ error 
(less than 25 percent) in precipitable water and only ±3.3 m error in the 500 to 920 mb 
thickness, using “locally” conditioned regression coefficients. The inadequacies of the 
“globally” conditioned retrievals are most dramatically displayed in retrieved thickness 
fields. 

This prelaunch study of VAS rctrievability gives promising results for pre-thunderstorm en- 
vironments. Profile errors are comparable to those claimed for pokir orbiting sounding instru- 
ments. Strong gradients are discernible in realistic partly cloudy conditions, although their magni- 
tudes are underestimated like those observed with operational sounders. Wliile not all of the com- 
plex factors in satellite data have been simulated in this study (such as multiple cloud structures, 
surface irregularities, reflected sunlight, and radiometric errors), these results still imply that VAS 
radiances will produce the best sounding when a human being classifies the scene, picks relatively 
clear areas for retrieval, and applies a “local” statistical database. The VAS simulated soundings 
show the same limitations in absolute accuracy and vertical resolution that are suffered by other 
passive infrared satellite instruments, especially in the lower tropospheric moisture retrievals. 

This weakness, even in these optimally simulated soundings, indicates that timely incorporation of 
conventional surface data will have to become part of VAS data processing to make the lower 
tropospheric retrievals more accurate. Nevertheless, VAS’ spatial and temporal resolution should 
permit us to locate and effectively monitor the relative changes within potentially unstable meso- 
scale structures, and thus significantly improve upon the synoptic scale, ground-based and polar 
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orbiting data networks. The analysis of real VAS data will determine if this potential for moni- 
toring mesoscale atmospheric changes can be realized with geosynchronous satellite sounding data. 
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TABLE 


Prelaunch specincations lor the 12 VAS channels, all of which have 15 km (mulir view) resolution. 


VAS 

CH 

# 
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CFNT 

n 

FILT 
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1 
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10 
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CO 2 
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.3 

38 
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691 

16 

temp 

CO 2 

0.4 

2.2 

.3 

10 

3 

14..^ 

699 

16 

temp 

CO, 

O 3 

1.8 

.3 

7 

4* 

14.0 

713 

20 

temp 

COj 

O 3 

1.2 

««v 

4 

S* 

13..1 

750 

20 

temp 

C 02 

HjO 

l.O 


4 

6 

4.5 

2209 

45 

icinp+cloud 

N,0 

sun 

1.6 

.1 

35 

1 * 

12.7 

787 

20 

moisuirc 

ILO 

CO, 

*1. 

LO 

.fc. 

4 

8 % 

11.2 

892 

140 

surface 


ILO+CO, 

0,1 

.2 

1 

()* 

7.2 

1370 

40 

moisture 

11,0 


3.4 

.4 

25 

lO* 

6.8 

1467 

150 

moisture 

11,0 


1.6 

.5 

-1 

II 

4.4 

2254 

40 

temp+cloud 

CO, 

sun 

6.7 

.3 

46 

12 

3.9 

2540 

140 

surface 


SU 11 +II 2 O 

0.8 

.1 

4 


*Available at 7 km (nadir view) resolmion. 
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TABLE 2 

Estimated sensitivity of the VAS channels for nadir view simulations of the U,S. Standard 
Atmosphere and a set of 32 “globally” distributed radiosondes. 


VAS 

CH 

# 

PEAK 

WEIGHT 

(mb) 

10%-90% 
RADIANCE 
LAYER (mb) 

NET SENSITIVITY 
CKI X) 

TEMP DEWPT SURF 

EARTH 

MEAN 

(“K) 

EARTH 

(±°K) 

“SIGNAL 

/NOISE 

CKI°K) 

1 

70 

4-120 

1.00 

-.00 

.00 

219 

10 

2 

2 

125 

4-225 

.99 

-.00 

.00 

219 

9 

4 

3 

200 

15-475 

.91 

-.02 

.00 

223 

7 

4 

4 

500 

40-900 

.73 

-.09 

.02 

238 

9 

7 

5 

920 

420-surf 

.47 

-.30 

.29 

260 

15 

15 

6 

850 

520-surf 

.43 

-.03 

.44 

257 

17 

11 

7 

1000 

720-surf 

,33 

-.37 

.63 

270 

19 

19 

8 

surf 

920-surf 

.17 

-.20 

.84 

272 

21 

210 

9 

600 

110-820 

.82 

-.70 

.02 

254 

10 

3 

10 

400 

240-620 

.94 

-.70 

.00 

246 

7 

4 

11 

300 

7-1000 

.80 

-.01 

,04 

237 

11 

2 

12 

surf 

surf-surf 

.08 

-.03 

.92 

273 

22 

27 
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TABLE 3 


Tile NSSL6 datasets used to “locally” condition regression matrices for simulated VAS soundings. 


FRAME TIME DATE 

1 2030 GMT 12 May 1976 


2030 GMT 21 May 1976 


3 1900 GMT 26 May 1976 


4 1900 GMT 12 June 1976 


5 2330 GMT 14 June 1976 


6 0100 GMT 18 June 1976 


DESCRIPTION OF 6 NSSL EVENTS 
A strong east-west squall line, associated with 
cyclogencsis in the Texas panhandle, moved to 
the south through the network. 

Scattered storm cells near the northwest corner 
of the network were associated with a stalled 
lecside trough west of Oklahoma. 

Areas of scattered activity developed over the 
network near the warm and cold fronts associ- 
ated with a low pressure system tliat moved 
over the area from the Texas panhandle. 

A few cells appeared in the northwest corner of 
the network just ahead of a slowly propagating 
dryline. 

A strong storm developed in response to a short 
wave propagating southeast over the network, 
initiated by the intersection of ,a frontal boun- 
dary with a pre-existing dryline, 

Strong storms developed in response to frontal 
and dryline interaction, as a weak short wave 
propagated eastward through the network. 
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TABLE 4 

The NSSL4 datasets used to test simulated VAS soundings. The four frames are drawn 
from the initial and developing stages of two well observed storms. The frame at 2030 
GMT 29 May 1976 is used to simulate VAS images. 


FRAME 

TIME 


DATE 


7 

2330 

GMT 

22 May 

1976 

8 

0100 

GMT 

23 May 

1976 


9H< 

2030 

GMT 

29 May 1976 

0 

2330 

GMT 

29 May 1976 


DESCRIPTION OF 2 NSSL EVENTS (4 frames) 
A strong upper level wave broke down a block- 
ing ridge and inversion to induce dynamics that 
resulted in a dryline frontal passage. A strong 
squall line developed along a moist tongue ahead 
of the front. 

An upper levei long wave over the Rocky Moun- 
tains supported cyclogenesis in the Texas-Okla- 
homa panhandles. A weak cold front combined 
with a wedge of hot, dry air behind a dryline 
and rapid moisture inflow from Gulf of Mexico, 
providing a classic example of rapid convective 
development directly over the NSSL network. 


*Cloud cover estimate is available from SMS 1 1 micron image. 
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TABLE 5 


Mean and sfaiulard deviations of prccipitable water and geopotential thiekness frotn the 
NSSL4 test set lor the “ground truth” profiles and residual errors in the “global” and 

“loeal” retrievals. 


CASE 


“Ciround 1'ruUi” 
“CJlobal” KesUluals 
“Local” Residuals 


RRECI STABLE WATER 
(gni enr“ , i 00-920 mb) 

2.49 ± 0.39 

-O.OI ± 0.39 

-O.IO * 0.I6 


GEOPOTENTIAL THICKNESS 
Oneters, 500-920 mb) 

4998. ±28.5 

-30.5 ± 6.9 

+ n I + IK 
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Figure 1. VAS weighting functions for 


Ir/d in P 


the seven temperature sounding channels, 
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Figure 3. SMS 1 1 micron image on 29 May 1976 at 2030 GMT with an overlay of the 21 x 1 1 
Barnes analysis grid at 20 km resolution for simulated VAS images and an imbedded 7x5 grid 

at 50 km resolution for simulated soundings. 
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Figua* 4b. Synthetic VAS briglitncss temperature images for a pre-storm environment at 2030 
GMT 29 May 1976, using the 20 km analysis grid and SMS cloud estimates. 
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RMS TEMPERATURE (K) 


Figunj RMS temperature residuals (simulated souudiiig »» “ground truth”), where VAS retrievals 
were statistically conditioned to expect vertical structure determined by either; (A) the “global” 

NOAA32 or (□) the “local” NSSL6 radiosonde sots. 
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Figure 7. RMS mixing ratio residuals, where VAS retrievals were statistically conditioned to expect 
vertical structure determined by either: (A) the “global” NOAA32 or (□) the “local” NSSL6 radio- 
sonde sets, 
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0.00 1.00 2.00 3.00 4.00 5.00 6.00 7.00 


RMS TEMPERATURE (K) 

Figure 8. RMS temperature residuals using: (□) all, (A) “wet” and (O) “dry” SFOVs from the 
“local” NSSL6 training set to test even more detailed statistical conditioning of the simulat- 
ed VAS soundings. 
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RMS MIXING RATIO (GM/KG) 


Figure 9. RMS mixing ratio residuals using: (□) all, (A) ‘Vet” and (O) “dry” SFOVs from 
the “local” NSSL6 training set to test even more detailed statistical conditioning of the 

simulated VAS soundings. 
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Figure 1 0. NWS surface analysis depicting pressure (mb) and 

GMT 30 May 1976. 
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s^OOOO GMT 30 MAY 1976 


